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SU(2) Gluodynamics and HP^ cr-model Embedding: 
Scaling, Topology and Confinement 
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We investigate recently proposed HP^ cr-model embedding method aimed to study the topology of 
SU(2) gauge fields. The HP^ based topological charge is shown to be fairly compatible with various 
known definitions. We study the corresponding topological susceptibility and estimate its value in 
the continuum limit. The geometrical clarity of HP^ approach allows to investigate non-perturbative 
aspects of SU(2) gauge theory on qualitatively new level. In particular, we obtain numerically precise 
estimation of gluon condensate and its leading quadratic correction. Furthermore, we present clear 
evidences that the string tension is to be associated with global (percolating) regions of sign-coherent 
topological charge. As a byproduct of our analysis we estimate the continuum value of quenched 
chiral condensate and the dimensionality of regions, which localize the lowest eigenmodes of overlap 
Dirac operator. 
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I. INTRODUCTION 

The topological aspects of gauge theories had always 
been the prime topic for the lattice community. How- 
ever, the real breakthrough here is due to the construc- 
tion of chirally symmetric overlap Dirac operator on the 
lattice jj, which via Atiyah-Singer index theorem 
applied in lattice settings allowed to investigate the 
topology of equilibrium vacuum fields. Since then a lot 
of remarkable results were obtained, among which the 
discovery of global (percolating) regions of sign-coherent 
topological charge 0, 0| is worth to be mentioned. Note 
that the overlap-based approach to the gauge fields topol- 
ogy is rather involved technically, however, the intricacy 
of the method is not the only problem here. What is 
more important for us is the factual absence of geomet- 
rically clean microscopic picture behind the topological 
fluctuations in the overlap approach, which is in a sharp 
contrast with what is known from seminal papers Q 
about the topology of Yang-Mills fields (see also (3|). 

Recently the alternative definition of the topological 
charge |^ has been put forward, which is based essen- 
tially on the same topological ideas as were used in |^ Q ■ 
The essence of the method consists in finding the nearest 
(in the configuration space) to given SU(2) gauge back- 
ground configuration of scalar fields representing non- 
linear cr-model with target space being the quaternionic 
projective space HP" (see Ref. for missing details). 
The relevance of HP" spaces is by no means accidental, 
the close connection between SU{2) Yang^ills theory 
and HP" cr- models is known for edges Generi- 
cally the idea is to realize the original SU (2) gauge holon- 
omy assigned to every lattice link as the non-Abelian ge- 
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ometrical phase factor associated with motion of scalar 
particle possessing non-trivial internal symmetry space. 
However, except for classical gauge backgrounds it is not 
possible to achieve this exactly, hence the appearance 
of the term 'nearest configuration' above. Unfortunately, 
this is the weakest aspect of HP" tr-model approach since 
the notion of nearest configuration could not, in fact, be 
conceived rigorously in the case of quantum fields and we 
have to rely mostly on the results of numerical simula- 
tions. However, we note that this is not a fatal restriction 
of the method, essentially the same problem exists, for 
instance, in overlap-based approach, where it was by no 
means evident a priori that chirally symmetric Dirac op- 
erator applied to thermalized fields is indeed local. Note 
that for equilibrium configurations the dependence upon 
the CT-model rank n was found @ to be trivial [s^ and 
for this reason we consider only the case n = 1 below. 

The crucial advantage of HP^ approach is its clean ge- 
ometrical meaning, which not only allows to capture the 
gauge fields topology, but also permits the microscopic 
investigation of the topological fluctuations - the luxury 
which was not available so far. The prime purpose of the 
present publication is to put this assertion on numeri- 
cally solid grounds. Note that some remarkable results 
we obtained already in Ref. |^ and we by no means cast- 
ing doubts on them. However, as far as the phenomeno- 
logically relevant quantities are concerned, the numerics 
should evidently be improved and we do this in present 
paper. In particular, after short theoretical introduc- 
tion and description of the numerical algorithms (sec- 
tionni) we present in section statistically convincing 
comparison of HP^ and overlap-based topological charges 
both measured on the same thermalized vacuum config- 
urations. Furthermore, in section Hvl the scaling of HP^ 
based topological susceptibility with lattice spacing and 
volume is discussed in details and its rather unambiguous 
extrapolation to the continuum limit is presented. The 
outcome of this calculations makes us confident that HP^ 
topological charge is fairly compatible with overlap-based 
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definition and is insensitive to lattice dislocations. More- 
over, keeping in mind the technical complexity and geo- 
metrical obscurity of the overlap construction, we dare to 
claim that the HP^ cr- model embedding method is indeed 
advantageous both geometrically and computationally. 

The geometrical explicitness of HP^ approach is cru- 
cial for the remaining part of our paper. Namely it allows 
us to investigate the topology related content of the orig- 
inal gauge fields in explicit manner by introducing the so 
called HP^ projection to be discussed in details in sec- 
tionEl Note that the term 'projection' here has nothing 
in common with its nowadays conventional meaning, nev- 
ertheless we use it mostly due to the historical reasons. 
The results we obtain this way are indeed remarkable, 
however, before mentioning them we would like note the 
following. In this paper we intentionally leave aside any 
serious attempts to interpret the experimental data. In 
fact, the theoretical interpretation could indeed be given, 
moreover, it overlaps parti ally with moderns trends in 
the literature [IS Hi E iH Q El- However, con- 
sideration of the theoretical implications would leads us 
to quite lengthy discussions and for that reason we de- 
cided to focus exclusively on the experimental side of the 
problem. Our theoretical interpretation of the results 
presented here is to be published elsewhere [l^ . 

Section is devoted to the investigation of the dy- 
namics of HP^ projected fields obtained from vacuum 
configurations at various lattice spacings. We start in sec- 
tion from the consideration of the curvature of pro- 
jected fields, which turns out to be extremely small. Fur- 
thermore, its spacing dependence does not contain any 
sign whatsoever of the perturbation theory contribution. 
Instead, the curvature vanishes in the limit a —^ 0, the 
leading spacing corrections being quartic and quadratic 
in a. In turn, these leading terms are to be interpreted 
as the gluon condensate and quadratic power correction 
to it, for both of which we obtain numerically precise 
estimations and compare them with available literature. 
Then in section FVBI the confining properties of HP^ pro- 
jected fields are investigated from various viewpoints. 
The conclusion is that these fields are indeed confining, 
moreover, the corresponding string tension accounts for 
the full SU (2) string tension in the continuum limit. The 
geometrical clarity of the HP^ approach allows us to in- 
vestigate in section IVB 21 the microscopic origin of the 
projected string tension. We present clear evidences that 
the string tension is to be associated with global (perco- 
lating) regions of sign-coherent topological charge men- 
tioned above. Finally, to conclude the study of HP^ pro- 
jected fields dynamics we consider in section FV CI the low 
lying spectrum of overlap Dirac operator in the projected 
background and show that both the spectral density and 
the localization properties of low lying modes are essen- 
tially the same as they are on the original fields. The 
quenched chiral condensate and the dimensionality of lo- 
calization regions as they are seen within HP^ projection 
are estimated. 

The results of section could only be convincing pro- 



vided that we make sure that the gauge fields compo- 
nents not surviving HP^ projection do not contribute to 
neither of non-perturbative observables discussed above. 
The problem is non-trivial due to the explicit gauge co- 
variance of HP^ projection, however we believe that it is 
addressed adequately in section IVTI where we show that 
this part of the original gauge fields could be defined 
unambiguously. Then it is the matter of straightforward 
calculation to show that the correspondingly constructed 
configurations are: 

i) topologically trivial, in a sense that the topologi- 
cal charge in both HP^ and overlap definitions vanishes 
identically on all configurations we have in our disposal; 

ii) deconfining, the string tension vanishes while the 
heavy quark potential is still compatible with Coulomb 
law; 

in) trivial from the fermionic viewpoint, the low lying 
spectrum of the Dirac operator disappears completely, 
we were unable to find even a single low lying eigenmode 
on all our configurations. 

These findings are the most stringent evidences that 
what had been left aside by HP^ projection corresponds 
to pure perturbation theory. At the same time, they 
show that the results and estimations of various non- 
perturbative observables, which were obtained in sec- 
tion are indeed reliable. Keeping in mind that the 
actual numbers are in complete agreement with the ex- 
isting literature, we dare to claim that the HP^ projec- 
tion is, in fact, the unique method allowing to achieve 
numerically most accurate results. 

Finally, let us mention that the approach developed in 
sections |V| EJ is parallel in spirit to what is known in 
the literature as 'vortex removal' [17j and 'monopole re- 
moval' ,lSj procedures. However, we stress that the HP^ 
projection is by no means similar to them both physically 
and technically. 



II. THEORETICAL BACKGROUND AND 
NUMERICAL ALGORITHMS 

In this section we briefly outline the recent approach 
aimed to investigate SU{2) gauge fields topology. This 
material is included to introduce the necessary notations 
and to make the paper self-contained (for details and 
further references see |5|). In particular, in this section we 
precisely describe the numerical methods utilized in our 
paper. Note that the present numerical implementation 
is slightly different from that of Q although we checked 
that the difference is completely inessential physics- wise. 

Quatcrnionic projective space HP^ can be viewed 
as the factor space HP^ = S"* = S'^ / S"^ (second 
Hopf fibering) and its lowest non-trivial homotopic 
group is 7r4(IIP^) ~ '^siS'^) ~ Z. The simplest ex- 
plicit parametrization of HP^ is provided by normalized 
quatcrnionic vectors 

k> = [go,9i]^, 9.eH, = g^g, 1 e H, (1) 
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where [■■■]'^ indicates transposition, H is the field of real 
quaternions, = gfcc S H, a = 0, ...,3 and Cq, are the 
quaternionic units 

eo = l, CiBj = ~6tj ~ EijkBk , hj,k^ 1,2,3 (2) 

with conjugation defined as usual 

qi^q°'ea, eo = eo , Ci ^ -a , i = l,2,3. (3) 

The states \q) describe 7-dimensional sphere, \q) & S"^ , 
while the HP^ space is the set of equivalence classes of | g ) 
with respect to the right multiplication by unit quater- 
nions (elements of SU(2) group) 



ID'", 



|2 — 



vv — 1 , u G H . 



(4) 



The gauge invariance I0J allows to introduce almost 
everywhere inhomogeneous coordinate uj = qiqQ^ = 
y°'ea e H on HP^ space 



(5) 



The alternative parametrization is provided by 5- 
dimensional unit vector n^, A = 0, ...,4 defined by 



(6) 



where 7^ are five Euclidean Dirac matrices {7", 7^} 
viewed as 2 x 2 quaternionic ones. One could show that 
© is equivalent to the standard stereographic projec- 
tion which relates n'^ with inhomogeneous coordinate 
u! = y°'ea above 



1 + r 



, a = 0, , 



1-r 



(7) 



It is clear that non- linear a-model with target space HP^ 
would be non-trivial provided that the base space is taken 
to be 4-sphere S'^. The corresponding topological charge 



Q 



(8^)2 



d X s 



£abcde ■ 



(8) 



is geometrically the sum of oriented infinitesimal volumes 
in the image of n'^{x) : S** — > HP^. 

As usual it is convenient to introduce auxiliary SU(2) 
gauge fields 



A^^-Al^-{q\d^\q), 



(9) 



transforming as A^ — > vA^v — vd^v under Q). Then the 
topological charge ijSJ) is expressible solely in terms of A^ 

Q = ^ y e^^'^^P Tr F,,Fxp , (10) 

= dfiAi, di/Afj^ -\- [Afj^, Au\ , 

being essentially equivalent to the familiar topological 
charge of the gauge fields ©. Note that Eq. (|10|l is by 



no means accidental, the deep connection between HP" 
CT-models and SU(2) Yang-Mills theory is known for edges 
(see, e.g., Refs. particular, the geometry 

of gauge fields could best be analyzed in the HP" CT- 
models context. In fact, all known instantonic solutions 
of Yang-Mills theory could be considered as induced by 
the topological configurations of suitable HP" CT-models. 

As is discussed in length in Ref. ^ , it is natural to ex- 
pect that the construction of nearest (in the configuration 
space) to the given gauge background HP" fields cap- 
tures accurately the gauge fields topology leaving aside 
the non-topological properties of the background. Note 
that for classical gauge fields the corresponding HP" ct- 
model reproducing Eq. © is known to be unique 0. 
However, the mathematical rigour is lost in case of 'hot' 
vacuum configurations and we could only hope to find 
the nearest in the sense of Eq. lO CT-model fields. The 
proposition of Ref. was to minimize the functional 



F{A,q) = {A, + {q\d,\q)f 



(11) 



for given A^ with respect to all possible configurations 
oi\qx)- On the lattice this equation translates into 



F{U,q) 



1 



— VSc 



\ Qx I qx+fi ) 

\{qx\qx+^)\ 



(12) 



where lattice gauge fields are denoted by Ux,^ € SU{2), 
V is the lattice volume and quaternionic scalar part is 
Sc[q] — {q + q)/2. The configuration | q^ ) which provides 
the minimum to F{U, q) is the best possible HP" fields 
approximation to the original gauge potentials 



U « C/"P = 



{ Qx I qx+^l ) 

\{qx I qx+,i )\ 



(13) 



while the minimum value of (|12|) is a natural measure of 
the approximation quality. Furthermore, it was shown 
in Ref. ^ that for equilibrium gauge configurations the 
dependence on the CT-model rank n is trivial since the 
quality of the approximation IjlSI) docs not depend on 
n. Therefore, it seems to be legitimate to consider solely 
the case n = 1 and only HP^ CT-model embedding is in- 
vestigated below. Note that usually the extremization 
tasks like (fl^ . (fT^ are expected to be plagued by Gri- 
bov copies problem. However, basing on the results of 
Ref. 5] we assert that this issue is most likely irrelevant. 
Additional indirect evidence is provided by the fact that 
we changed slightly the extremization algorithm (see be- 
low) compared to that of Ref. however, all physically 
relevant results turned out to be the same. As far as 
the actual numerical algorithms used in this paper are 
concerned, they could be summarized as follows. 

• HP^ CT-model embedding. 

The functional (|12|l was minimized sequentially at each 
lattice site, the non-linearity was taken into account by 
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performing 3 iterations 





a,fm 


Lt 


Ls 









(14) 



at each point x and then sweeping number of times 
through aU the lattice until the algorithm converges. The 
proportionality factor, which is omitted in (|14|l . is chosen 
to ensure the normalization condition ( qi^ \ qi^ ) = 1 at 
each iteration. The stopping criterion was that the dis- 
tance between old and new | q ) values at all lattice sites 



max (1 



Note that the value of r.h.s. is the result of trade-off 
between the computational demands and the needed nu- 
merical accuracy (see below). It turns out that Eq. ifT^ 
provides the accuracy of miug F(U, q) calculation of or- 
der 10^^ which is the same as was used in 13. As far 
as the Gribov copies problem is concerned, we always 
considered 10 random initial distributions | qx ) and then 
selected the best minimum found. 

• Global topological charge and topological density. 

The topological charge density in terms of the HP^ a- 
model fields is given by the oriented 4-volumes of spher- 
ical tetrahedra T embedded into HP^ = (image of 
n"^ : iS*^ HP^). Unfortunately, the only way to esti- 
mate these volumes is to use the Monte Carlo technique 
which however forbids the exact volume evaluation. In 
this paper we used precisely the same algorithm of topo- 
logical chM|ge density calculation which was discussed in 
length in In particular, the embedded HP^ cr-model 
assigns 5 unit five-dimensional vectors nf i = 0, ...,4 to 
every simplex T of physical space triangulation and the 
topological charge density in simplex T is given by 

q{T) = ^ sign (det,^ [n^)]) • V{T) , (16) 

where V{T) is the volume of spherical tetrahedron T and 
normalization factor 8tt'^/3 is the total volume of S*. The 
corresponding topological charge 



Q float = q{T) 



(17) 



T 



is not integer valued due to finite accuracy of Monte Carlo 
estimation of V{T). However, for high enough accuracy 
of calculation at each simplex the non-integer valuedness 
of the topological charge is in fact irrelevant and we could 
safely identify 



Q = [Qf 



loat I 



(18) 



where [x] denotes the nearest to x integer number. The 
validity of the identification H18|) was discussed in details 



Fixed Volume 



2.3493 


0.1397(15) 


10 


10 


3.8(2) 


300 
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2.3772 


0.1284(15) 
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10 


3.8(2) 


300 
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2.3877 


0.1242(15) 


16 


10 


3.8(2) 


300 




2.4071 


0.1164(15) 


12 


12 


3.8(2) 


200 




2.4180 


0.1120(15) 


14 


12 


3.8(2) 


100 




2.4273 


0.1083(15) 


16 


12 


3.8(2) 


250 




2.4500 


0.0996(22) 


14 


14 


3.8(2) 


200 




2.5000 


0.0854(4) 


18 


16 


3.92(7) 


200 





Fixed Spacing 





2.4180 


0.1120(15) 


14 


14 


6.0(3) 


200 




(15) 


2.4180 


0.1120(15) 


16 


16 


10.3(6) 


200 






2.4180 


0.1120(15) 


18 


18 


16.5(9) 


200 





2.4000 


0.1193(9) 


16 


16 


13.3(4) 


198 




2.4750 


0.0913(6) 


16 


16 


4.6(1) 


380 


40 


2.6000 


0.0601(3) 


28 


28 


8.0(2) 


50 





TABLE I: Simulation parameters. 



in 13 . Here we note that there exist an additional cross- 
check of Eq. H18|l based on the observation that in order 
to get the global topological charge it is not necessary to 
calculate its density in the bulk. It is sufficient to count 
(with sign) how many times a particular point tiq € HP^ 
is covered by the image of : ^ HP^ 



E 



signdetiA [nf^^] 




if no e T 
otherwise 



(19) 



Evidently, Eq. (|19|) is computationally superior to H18|) 
and we used it in present paper to calculate the topo- 
logical charge. Moreover, the comparison of (|18l) and 
(|19|l provides the stringent calibration of Monte Carlo 
topological charge density evaluation It goes with- 
out saying that we indeed checked that both definitions 
(|18|l . (|19|l give exactly the same topological charge on all 
available configurations. Finally, we remark that slicing 
of hypercubical lattice into simplices was done according 
to 19]. 

• Simulation parameters. 

Our numerical measurements were performed on 14 
sets (Table IJ) of statistically independent SU(2) gauge 
configurations generated with standard Wilson action. 
The data sets are subdivided naturally into three groups 
indicated in Tabled The fixed volume configurations are 
exactly the same as ones used in Refs. 13, Tl^ in studying 
low lying Dirac eigenmodes localization properties. This 
allows us to make statistically significant comparison of 
HP^ cr-model approach and the overlap-based topological 
charge definition (section llll|) . Gauge configurations at 
fixed spacing were generated to investigate the finite vol- 
ume dependence of HP^-based topological susceptibility 
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f section Hvjl . In addition to the gauge configurations at 
P = 2.400, 2.475 which were analyzed already in Ref. 0, 
we also include the data set a.t f3 — 2.600 with finest 
lattice spacing we have so far. 

The last column in Table represents the number of 
configurations on which we calculated the bulk topolog- 
ical charge density using the above described algorithm. 
Note that it is slightly different from that of Ref. 
and this explains the relatively small number of analyzed 
configurations. We stress that the results obtained with 
modified method remain essentially the same. Neverthe- 
less, for honesty reasons the old data is not included in 
the present study. The relatively small number of config- 
urations on which we calculated the bulk topological den- 
sity distribution does not spoil the statistical significance 
of our results. Indeed, the majority of observables to be 
discussed below depend only upon the global topological 
charge obtained via Eq. (|19|l and upon the corresponding 
HP^ cr-model. In turn the HP^ fields were constructed 
for all configurations listed in Table ^ thus ensuring the 
statistical significance of our results. 

The lattice spacing values quoted in Tabled were par- 
tially taken from Refs. |2(1 |2l| and fixed by the physical 
value of SU(2) string tension y/a = 440 MeV. Note that 
not all /3-values listed in Table^could be found in the lit- 
erature. In this case the lattice spacings and correspond- 
ing rather conservative error estimates were obtained via 
interpolation in between the points quoted in [20l I2H . 



III. Hpi-BASED TOPOLOGICAL CHARGE OF 
SEMICLASSICAL AND VACUUM FIELDS 

The validity of HP ^ -based topological charge definition 
could only be convincing provided that we confront it 
with other known topological charge constructions. The 
preliminary study of this problem was already under- 
taken in Ref. however, only correlation of local topo- 
logical densities in HP^ and overlap-based approaches on 
relatively small statistics was considered in that paper. 
Here we fill this gap and compare the distribution of the 
HP^ global topological charge Qhp^ , Eq. II19II . with that 
in field-theoretical 22] and overlap-based jll HI construc- 
tions. 

Before presenting the results of our measurements let 
us discuss what could be the proper quantitative charac- 
teristic of the correlation between various definitions of 
lattice topological charge. The problem is that there is no 
commonly accepted viewpoint in the literature, usually 
people just consider the cumulative distribution of the 
topological charge in pair of definitions. For instance, in 
Ref. '23| the topological charge distribution histogram in 
the plane .56] {Qdover, Qoveriap) was presented, however, 
no quantitative measure of correlation between Qdover 
and Qoveriap was given. Below we discuss two observ- 
ables which quantify the notion of correlation between 
the topological charges Qi and Qj obtained via some ab- 
stract constructions "i" and "j" correspondingly. The 



first one could be obtained from the above mentioned 
cumulative distribution in {Qi,Qj) plane by noting that 
the complete equivalence of "i" and "j" constructions 
would give Qi — Qj identically and therefore the en- 
tire plot would collapse into Qi — Qj line. On the other 
hand, any disagreement of "i" and "j" approaches would 
make the cumulative histogram broader and therefore it 
is natural to consider the ratio of the widths of cumula- 
tive distribution projected onto the lines Qi — ~Qj and 
Qi — Qj I for which the limiting values are (exact equiv- 
alence of "i" and "j" constructions) and ~ 1 (no correla- 
tion of Qi an d Qj ) . Moreover, it is similar to what had 
been done in |23| thus allowing to compare our results 
with the existing literature. However, it is clear that the 
widths ratio is mostly graphical characteristic since the 
distribution along Qi — —Qj line is not guaranteed to be 
Gaussian-like, hence the corresponding width is ambigu- 
ous in general. To ameliorate this we propose to consider 
the correlator 



{Q^ Ql) 



(20) 



where averages are to be taken on the same set of configu- 
rations. Note that the normalization in (|20|) is such that 
rfj K, 1 for essentially equivalent "i" and "j" construc- 
tions, while rfj = in the case of complete decorrelation 
of Qi and Qj. 

With all these preliminaries, let us compare our ap- 
proach with field-theoretical method, which in turn could 
be applied reliably only on semiclassical (cooled) config- 
urations. Thus its comparison with Qup^ is rather tech- 
nical issue related mostly to the validity of Qup^ con- 
struction on semiclassical background. To this end we 
generated 300 cooled configurations (not included into 
TablePi initially thermahzed at /3 = 2.400 on 16^ lattice. 
The cooling algorithm em ploy ed in our study is anal- 
ogous to one described in j2J]. The original fields were 
cooled until the action stabilizes indicating that semiclas- 
sical regime had been reached. At this point we applied 
the algorithm of HP"'^ cr-model embedding thus obtaining 
Qnpi and simultaneously measured the field-theoretical 
topological charge Qdover defined by 

Qdover — Qdover{-^^ : (^-i) 

X 

Qdoveri.-^^ — 327f^^/^i^Ap Tr [ ^^^^ (a^) _F\p (x) ] , 

where F,,,^ is the 0{a'^) improved lattice field-strength 
tensor [23, H^- The resulting cumulative distribution 
of Qhri Eind Qdover is presented on Fig. ^ As is evi- 
dent from that figure both definitions are fairly compat- 
ible with each other making us confident that Qjjpi ^^'^ 
Qdover identify essentially the same topology on all con- 
sidered data sets. The disagreement between Qnpi and 
Qdover is Only secn in a few points which, however, is to 
be expected (see Ref. |23] for discussion of this issue in 
case of field-theoretical and overlap constructions). For 
this reason the above discussed widths ratio is certainly 
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FIG. 1: Cumulative distribution of Qhpi i Eq. 1191 and 
Qciover, Eq. 12111 . on 300 cooled configurations initially ther- 
malized at /9 = 2.400 on 16^ lattice. 
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FIG. 2: Cumulative distribution of Qhpi i Eq. 11911 and 
Qoveriap, Eq. 124L ou 200 fixed volume configurations at 
P = 2.500, Table 13 



very small and is comparable with zero within the sta- 
tistical errors. As far as the correlator (|20(l is concerned, 
its numerical value is 

T^'^^T'' = 0.94(1) (22) 

and indeed is very close to unity. 

Turn now to the most interesting comparison of HP^ 
topological charge with becoming standard now construc- 
tion based on Atiyah-Singer index theorem applied to chi- 
rally symmetric overlap Dirac operator ^] . Explicitly 
the overlap operator is given by 

D=P(l + ^£=), A^Dw-^, (23) 

where A is the Wilson Dirac operator with negative mass 
term and we used the optimal value 1.4 of p param- 
eter. Anti-periodic (periodic) boundary conditions in 
time (space) directions were employed. To compute the 
sign function sign(A) — A/vAA^ = 75sign(iJ), where 
H = 75A is hermitian Wilson Dirac operator, we used 
the minmax polynomial approximation (211] . In order to 
improve the accuracy and performance about 50 lowest 
eigenmodes of H were projected out. Note that the eigen- 
values of D are distributed on the circle of radius p cen- 
tered at (p, 0) in the complex plane. Below we will need 
to relate them with continuous eigenvalues of the Dirac 
operator and therefore the circle was stereographically 
projected onto the imaginary axis p3.l2^. 

It is well known that the knowledge of the spectrum of 
(|23|l allows to define the topological charge via Atiyah- 
Singer index theorem 0| applied in lattice settings 3] 

Q overlap — ^-\- : (^4) 

where n+ (ji-) is the number of positive (negative) chi- 
rality zero eigenmodes. In order to compare Qnpi and 
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FIG. 3: Projection of the cumulative distribution of Qnpi ^-^d 
Qoveriaj] sX (3 = 2.500 outo the fines Qhpi = ^Qoveriap, aloug 
which Qproj \/2 is the corresponding integer valued coordi- 
nate. Note that both histograms are normalized separately. 



Qoverlap wc Considered both the cumulative distribution 
and the quantity 1)20(1 on 200+200 fixed volume configura- 
tions at /3 = 2.450 and P = 2.500 (see TablePl. Note that 
since the width of Q distribution is generically spacing de- 
pendent it would be unwise to mix Qnpi , Qoverlap cumu- 
lative distributions at different /3-values and we present 
on Fig. 121 only (3 = 2.500 results. Evidently both defi- 
nitions are strongly correlated with each other although 
the distribution is notably broader than that on Fig. 2] 
The relative broadness of Qnpi i Qoverlap distribution is 
not surprising by itself, clearly it must be broader than 
that on cooled configurations. What is relevant here is its 
projection onto the lines Q^pi = ^Qoverlap and the ratio 
of corresponding widths. The projected histograms are 
presented on Fig. |31 from which it is clear that they both 
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are Gaussian-like, the widths ratio being about 1/3-^1/4. 
To get the feehng of the numbers involved it is instruc- 
tive to remind the results of Ref . jl^l , where the cumula- 
tive distribution of Qdover and Qoveriap was considered. 
Note that the relevant plot of that paper refers to cooled 
configurations and its direct comparison with Figs. El 
would be unfair since for cooled configurations the distri- 
bution is generically expected to be even narrower. How- 
ever, it is pleasant to notice that the ratio of the corre- 
sponding widths in Qdover, Qnveriav distribution, which 
could be extracted from Ref. [23, is larger than that of 
QhpIj Qoveriap and is of Order 1/3-^ 1/2. Thus we con- 
clude that Qnpi J Qoveriap Correlation seems to be even 
stronger than it is in the case of Qdover and Qoveriap 
charges. Moreover, the strong correlation of Qhr^ and 
Qoveriap is also Confirmed by the quantity (|20|l and what 
is probably more important is that it is rising with di- 
minishing spacing 

»7hpT'"''(/3 = 2.450) = 0.71(1), 
^o^^eHap^^ ^ 2.500) = 0.77(1). 

Note however that unlike the previous case of widths ra- 
tio the numerical values of rj'^pi^'^^ cannot be compared 
with the existing literature. 

To summarize, we found clear and statistically signif- 
icant evidences that the construction of the topological 
charge based on HP^ cr- model embedding is fairly com- 
patible with field-theoretical definition in case of semi- 
classical configurations and is strongly correlated with 
overlap-based approach in the case of equilibrium vac- 
uum fields. In the next section we put the HP^ method 
under other tests, namely, consider the scaling properties 
of topological susceptibility defined via HP^ cr-model em- 
bedding. 
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FIG. 4: The topological susceptibility as a function of lattice 
volume both being expressed in physical units (fixed spacing 
configurations at f3 = 2.4180, Tabled. 
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IV. SCALING OF TOPOLOGICAL 
SUSCEPTIBILITY 

The basic topology related observable in pure Yang- 
Mills theory is the topological susceptibility which is for- 
mally defined as 

X = J d^x{q{0)q{x)) . (25) 

Here q{x) is the topological charge density which is equal 
to qdover (x), Eq. H21(l . in the naive continuum limit. The 
topological susceptibility, being the quantity of prime 
phenomenological importance, is, in fact, not well de- 
fined within the definition H25|l and is written usually as 

X = Ihn {Q^)/V (26) 

in the context of lattice regularization, where a and V 
are the lattice spacing and lattice volume respectively. 



FIG. 5: Topological charge distribution obtained at /3 = 2.475 
on 16* lattices (380 configurations). Curve represents the best 
fit to Eq. l|?7jl . 



The purpose of this section is to study the limit l|26|l and 
obtain the estimate of the topological susceptibility in 
the continuum using the HP^ definition of the topological 
charge Q = Q^pi ■ 

Due to the absence of massless excitations in pure 
glue theory the topological susceptibility is expected to 
reach its infinite volume limit exponentially fast with 
increasing lattice size. The characteristic length of ex- 
ponential fall off is dictated by lightest glueball mass 
(1.5 GeV)~^ « 0.15 fm and therefore we expect that 
finite- volume effects are negligible for our lattices. To 
confirm this assertion we measured the topological sus- 
ceptibility on fixed spacing configurations (Table at 
P = 2.4180 which cover rather wide range of lattice vol- 
umes. The result is presented on Fig.^and shows clearly 
that indeed finite size corrections are much smaller than 
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FIG. 6: Scaling of the HP^-based topological susceptibility 
with lattice spacing; line is the best fit to Eq. 12811 . 



the statistical errors. 

An additional cross-check of the negligibility of finite- 
size corrections conies from the observation [23, that 
the probability distribution of the topological charge in 
large yolume limit is to be Gaussian 



Pc, 



1 



-Q^/(2{Q^)) 



(27) 



We haye checked that the distribution Pq is indeed as 
expected for all our data sets. A particular illustration 
is provided by the data taken aX (3 — 2.475 on IG'* lattice 
where we have the largest statistics (Fig. O. Moreover, 
we were unable to find any statistically significant devi- 
ations of Pq from 1)27(1 on all our data sets. 

As far as the lattice spacing dependence of the suscep- 
tibility is concerned, it is not a priori evident that x{^) is 
not plagued by power divergences contrary to the case of 
overlap-based definition. Unfortunately, we're still lack- 
ing the theoretical arguments which could highlight the 
issue and should rely, in fact, only on the results of the 
numerical simulations. Note however that at least non- 
universal 0{a?) discretization effects are expected gener- 
ically. Therefore let us consider the topological suscepti- 
bility as a function of squared lattice spacing, Fig. |S1 As 
is clear from the figure the dependence x(^^) is totally 
compatible with linear one 



(28) 



and does not show any pathology in the limit of small 
lattice spacing. Of course, the dependence (|28() is not 
guaranteed to persist in the academical limit a —^ 0, but 
our data strongly disfavors this possibility. The conclu- 
sion is that the topological susceptibility in the HP^ ap- 
proach is not plagued by power divergences. Moreover, 
the best fit according to Eq. H28|l leads to 



which is totally compatible with the conventional value 
of topological susceptibility in the continuum limit (see, 
e.g., Ref. [23| and references therein). Note that in 

Ref. [3 we got xlitiap = 225(3) MeV on essentially 
the same configurations using overlap-based definition of 
the topological charge. We believe that this difference 
is completely inessential and should be attributed to the 
finite width of cumulative Qnpi , Q overlap distribution 
along the line Qep^ = ^Qoveriap, which we discussed in 
previous section. 

To summarize, our high statistics computation of the 
topological susceptibility, which was always considered as 
a testbed for all topological charge constructions, shows 
that the HP^ method is most likely to be free of any 
lattice related pathology. In particular, we do see al- 
most perfect scaling of x within our approach and it does 
not reveal any sign of necessity of multiplicative renor- 
malization. This provides us a stringent evidence that 
HP^ topological charge is not sensitive to lattice disloca- 
tions. As far as the physical results are concerned, the 
HP^ topological susceptibility being extrapolated to the 
continuum limit is perfectly consistent with the results 
of other investigations, most notably with overlap-based 
approach. Keeping in mind all subtleties and technical- 
ities involved in the construction of Qoveriap we could 
dare to claim that the HP^ cr-model embedding method 
is comparable with overlap based definition and is, in 
fact, much more advantageous computationally. More- 
over, it is not only the tool to investigate the topological 
aspects of equilibrium vacuum gauge fields. As we show 
below it naturally allows to study the role of topology 
related fiuctuations in the non-pcrturbative dynamics of 
SU(2) Yang-Mills theory. 



V. cr-MODEL INDUCED GAUGE FIELDS 

The essence of HP^ cr-model embedding method is the 
assignment of quaternionic valued scalar fields { | ) } to 
the original SU(2) gauge fields configuration. The crucial 
question here is the uniqueness of this assignment and, 
unfortunately, in the case of equilibrium gauge fields it 
seems to be impossible to investigate this issue analyti- 
cally. However, the results of Ref. and the ones pre- 
sented above suggest that non-uniqueness is most likely 
to be physically irrelevant. Therefore let us assume that 
the association of {l^^)} with a particular gauge con- 
figuration is indeed unique. It is crucial that the corre- 
sponding HP^ projection 



u. 



HP" 



( ix I qx+f, ) 
l( I qx+^i ) I 



(30) 



X 



1/4 



216(4) MeV, 



(29) 



is gauge covariant since under gauge rotations both Ux.^i 
and transform exactly in the same way. In this re- 

spect the HP^ projection is radically distinct from what 
is usually called 'projection' in the literature. The very 
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appearance of Eq. (|30|l was motivated only by considera- 
tion of the gauge fields topology and the word 'projection' 
here is just the unfortunate terminology. However, this 
term is in common use and we adopted it also. 

Eq. (|3U|) naturally allows to consider the properties of 

{U^^ } fields viewed as usual SU(2) matrices assigned 
to lattice links. The purpose of this section is to investi- 
gate the dynamics of HP^ induced gauge fields from var- 
ious viewpoints. In section IVAI we consider the simplest 
observable, the gauge curvature of induced potentials, 
and investigate its dependence upon the lattice spacing. 
Then the crucial question of confinement in the projected 
fields is studied in section lV Bl Finally in section lV ("I the 
properties of {U^^ } configurations are investigated with 
chirally symmetric overlap Dirac operator. Since this sec- 
tion is rather lengthy and presents quite remarkable re- 
sults, we conclude it with brief summary fsection lVPp . 

A. Gluon Condensate and Quadratic Corrections 

The immediate question about HP^ induced gauge 
fields is the corresponding curvature. On the lattice this 
amount to the consideration of the trace of the plaquette 
matrix 1 /2 Tr [/^^ constructed as usual from f/^J^ and 
it turns out that HP^ induced curvature is astonishingly 
smah, (l/2Tr[/^pi^ > o.97, on all configurations we 
have in our disposal. To illustrate this we note that the 
typical distribution of 1 — 1/2 Tr U^^ measured on HP^ 
projected configurations is radically different from that 
on the original fields and is well described by power law, 
which should be compared with usual almost exponen- 
tial tail in 1 — l/2TrJ7p distribution. It is amusing that 

1 — 1/2 TrC/^^ distribution looks similar to the lumps 
volume distribution reported in and most probably 
they are indeed closely related. In fact, it is a clear ev- 
idence that the dynamics of HP^ projected fields is to- 
tally different from that of original ones although it is still 
completely determined by the original Wilson action. 

In principle it is an interesting question how the above 
power law dependence changes with lattice spacing, but 
we'll not directly consider it. Instead let us focus on the 
spacing dependence of averaged plaquette {l/2TiU^^ ) 
for which rather accurate results could easily be ob- 
tained. We have measured the HP"'^ induced mean pla- 
quette (l/2TrC/^P ) at various spacings quoted in Ta- 
ble^ The result is presented on Fig.[7|and indeed is quite 
remarkable. The spacing dependence of (l/2Tr?7^^ ) is 
totally distinct from what could be expected for unpro- 
jected fields and as a matter of fact is astonishingly well 
described by simple polynomial 

1 - (l/2TrC/jfP') = asa^ + a^a'^ (31) 

throughout the whole range of spacings considered. The 
justify the assertion we plot on Fig. [7| the depen- 
dence of a-2-(l- (1/2 TrC/^P )) upon the squared lattice 




I ' ' ' ' ' 1 

0.1 0.2 0.3 0.4 0.5 0.6 

FIG. 7: Scaling of averaged plaquette P = (1/2 TrC/^^') 
constructed from HP^ projected gauge fields. Line is the best 
fit according to Eq. 

spacing and it is apparent that this dependence is indeed 
linear. In turn the best fit gives /n.d.f. = 1.2 and the 
optimal values of parameters are 

a2 = [61(3) MeV]2 , (32) 
a4 = 0.0271(10) GeV^. 

It is quite remarkable that the dependence H31() shown 
on Fig. 13 by solid line includes only positive powers 
of lattice spacing. This means, in particular, that 
(1/2 Tr ) contains no trace whatsoever of the pertur- 
bation theory contribution and is vanishing in the limit 
a ^ 0. Thus 1 - (l/2TrJ7^P') could be considered as 
the definition of non-perturbative contribution present in 
{as G^/tt) full, where we have used the conventional con- 
tinuum notations and subscript indicates that averaging 
is done in the full theory including perturbative series. 
We are in haste to add, however, that at the time being 
this definition looks completely ad hoc. Indeed, we can- 
not identify (l/2TrC/HP ^ ^j^g ^^ly 

non-perturbative 

part of the {as G^/vr) fuii until we analyze the content of 
original gauge fields not surviving HP^ projection H30|l . 
We will do this in section IVll and right now let us agree 
to treat (1/2 Tri/^^ ) as the genuine non-perturbative 
quantity. Then it is straightforward to relate the coeffi- 
cient a4, Eq. (|31|) . with gluon condensate (a^ G^/tt) in- 
troduced first in 31]. Indeed, discarding temporally the 
O(a^) term in (|31|) wc have 

1- (l/2TrC/HP') = aia\l3)^ (33) 

= a\(3)^{asG^/^), 

thus obtaining the following estimation of the gluon con- 
densate 

{asG^ln) = 0.066(2) GeV* . (34) 
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As far as the comparison with existing hterature is con- 
cerned, let us mention that the gluon condensate being 
the quantity of prime phenomenological importance is 
most frequently discussed within the SU (3) pure gauge 
theory. Here the most recent lattice measurements give 
{as /n)su{3) ~ 0.04 GeV albeit with large uncertain- 
ties (see, e.g., Ref. and references therein). This 
value should be compared with phenomenological one 
{asC^Msuis) « 0.012 GeV coming from SVZ sum- 
rules jsj. We see that even in the most phenomeno- 
logically relevant case of SU{3) gauge group the estima- 
tions of the gluon condensate vary widely and are, in fact, 
fairly compatible with H34() . If we turn now to the case 
of SU{2) gauge group, the results which we were able to 
find in the literature are again in accordance with our 
estimation and vary from (a., /tt) w 0.02 GeV"*, 
[13, to (a^GVTr) w 0.15 GeV, Hj (for review see, 
e.g., Refs. Isll and references therein). Therefore, we are 
confident that the value of the gluon condensate com- 
ing from HP^ cr-model embedding method agrees rather 
nicely with what is known in both SU(2) and SU(3) cases. 
Moreover, in view of the well known tremendous uncer- 
tainties of usual approaches we conclude that Eq. H34|) 
provides the most accurate estimation of the gluon con- 
densate in SU{2) gauge theory available so far. Note 
that the systematic errors involved in (|34|) are coming 
only from HP^ projection H30I) . However, as we argue in 
section IVll the systematic uncertainties are most likely 
to be vanishing thus justifying the estimate H34|l . 

Turn now to the quadratic correction term present 
in (|31|l . Generically it corresponds to the dimension 
2 condensate (A^jn) which was introduced, in fact, 
long ago j36j and became the subject of active devel- 
opment recently 0, [s^l (see, e.g., Refs. ^3 and refer- 
ences therein). Note however that the theoretical sta- 
tus of quadratic corrections in general and, in particu- 
lar, the status O(a^) correction to the gluon condensate 
is uncertain at present. Indeed, the quadratic correc- 
tion to {asC^/ir) was clearly seen in Ref. [s^l, however, 
the corresponding coefficient decreases with increasing 
number of subtracted perturbative loops |38j| and it is 
by no means evident whether 0{a?) term has the non- 
perturbative origin or comes from the higher orders of 
perturbation theory (see Refs. for clear and con- 

cise discussion). Generically our data seem to dis- 
favor the perturbative origin of the quadratic term (see 
also section IVl|) . However, as was argued in Introduction 
we are not in the position to interpret theoretically this 
problem and would like to postpone the corresponding 
analyzes Il6j . 

As far as the actual numbers are concerned, it is re- 
markable that the value of a2 coefficient, Eq. (I32|l . is un- 
expectedly small compared to the natural scale Aq^^, 
but nevertheless almost fits into the established bounds 
quoted in Ref. [s^l- Note that although the data for 
SU{?>) gauge group was presented in pJl: it seems for us 
that the parametric smallness of the quadratic correction 
term is quite generic and should also be valid in SU (2) 



case. Therefore, we are confident that the estimate (|32|) 
is not in contradiction with the most recent literature. 

To summarize, the measurements of the curvature of 
HP^ projected gauge fields performed in a wide range of 
lattice spacings reveal that these potentials are extremely 
weak compared to that of the original Yang-Mills theory. 
Moreover, the lattice spacing dependence of averaged 
HP^ plaquette is astonishingly well described by simple 
power law, Eq. (|^ . with only 0{a!^) and O(a^) terms 
present and with no sign whatsoever of the perturbation 
theory contribution. We argued that these terms are to 
interpreted as the non-perturbative gluon condensate and 
unusual quadratic power correction to {asC^/ir). The 
factual absence of perturbative tail in (1/2 Trt/^^ ) al- 
lows us to obtain numerically precise estimation of the 
gluon condensate, which agrees nicely with the existing 
literature. On the other hand, the measured magnitude 
of the quadratic term allows to investigate the issue of 
unusual power corrections in Yang-Mills theory on qual- 
itatively new level. As far as the systematic bias intro- 
duced by HP^ projection (|30|l is concerned, we postpone 
its discussion until section IVTI and only note here that it 
is likely to be inessential. 



B. Confinement in Induced Gauge Fields 

In this section we consider the crucial question of con- 
finement in HP^ projected gauge fields. In section IVB II 
we discuss the Wilson loops confinement criterium and 
the scaling properties of HP^ projected string tension 
cr^^ . Then we try to identify the objects which seem to 
be directly related to the non-vanishing cr^^ . Finally, 
in section lVB3l the Polyakov lines correlation function is 
considered. 



1. String Tension from Wilson loops 

In order to investigate the confinement properties of 
HP^ projected gauge fields H30|l we have measured the 
planar TxR Wilson loops using various data sets listed in 
TableQ] As might be expected the weakness of projected 
fields allows to obtain numerically rather accurate re- 
sults even without the conventional spatial smearing |39l | 
and hypercubic blocking j40l | of temporal links. Indeed, 
we have checked that the utilization of both these tech- 
niques does not change our results in any essential way. 
However, both these tricks are in fact mandatory to ex- 
tract the string tension of the original gauge fields to be 
compared to that on {U^^ } configurations. In order 
to make the algorithms, applied to the original and HP^ 
projected fields, coherent the smearing and hypercubic 
blocking was done in both cases. From {W{T, R)) the 
heavy quark potential V^^_(R) was extracted in stan- 
dard way (see, e.g. Ref. |4j| for details). In order to 
obtain the HP^ projected string tension the poten- 
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FIG. 8: —ln{W{T, R)) as a function of T at various R mea- 
sured at /3 = 2.600 on 28* lattice. 




FIG. 9: Heavy quark potential V^^ (R) extracted from Wil- 
son loops measured at /3 = 2.600 on 28* lattice. Solid curve 
is the fit to Eq. 



tial was fitted to 



(R) = const + - 
R 



(35) 



To convince the reader that despite of the U^^^ fields 
weakness the potential (i?) is indeed linearly rising 
we show on Fig. |S| the behavior of — ln(M^(r, R)) as a 
function of T at various R measured at /3 — 2.600 on 28'' 
lattice. The corresponding potential V^^ (R) is depicted 
on Fig. El Note that the errors bars here are smaller than 
the size of symbols. It is apparent that the potential 
has positive second derivative for i? < 4 which probably 
indicates the reflection positivity violation. However, we 
don't think that this is a real problem. Indeed, the HP^ 
projected gauge fields do not contain the perturbative 
contribution (see above) and already for this reason are 
not obliged to fulfill the usual requirements of reflection 
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FIG. 10: Scaling of HP^ projected string tension cr^^ com- 
pared to that of full SU (2) string tension a^'^^^^ . Lines are the 
best fits according to Eq. Il.'i6|l in the whole range of available 
spacings and in the region a < 0.12 fm. 



positivity. Additional argument comes from the intrinsic 
non-locality of HP^ projection. We could only hope to 
recover the refection positivity in the large distance limit 
and it is pleasant to note that indeed the potential does 
not look pathological for i? > 3 (i? > 0.2 fm in physical 
units). 

As a matter of fact these plots are typical for all our 
data sets. The only difference is that at smaller /3-values 
we don't see any sign of the reflection positivity viola- 
tion at small R, the potential is everywhere linear within 
the errors bars. However, the projected string tension 
depends non-trivially on the lattice spacing. In order to 
investigate cr^^ scaling properties we calculated the full 
SU{2) string tension a^^^"^^ on our configurations and 

considered the ratio of va^^ and va^^-^ . This trick 
allows us not only to cross-check the physical scale of our 
fixed volume configurations, but also to avoid the finite 
volume systematic uncertainties. The resulting plot of 
the ratio [cr^^ l(jSU{2)Y/'i with rather pessimistic error 
estimates as a function of lattice spacing is presented on 
Fig. ^1 It is clear that the scaling of string tensions ratio 
is fairly compatible with linear a-dependence 



Co 



ci • a . 



(36) 



Indeed, on general grounds we could exclude the possibil- 
ity of negative powers of a in H36|l while our data is prac- 
tically insensitive to the inclusion of O(a^) terms. Appar- 
ently three points corresponding to the largest spacings 
systematically deviate from (|36|l which can be attributed 
the closeness of the crossover transition. Therefore we 
could try to fit the data to Eq. (|36|l either in the re- 
gion a < 0.12 fm or in the whole range of available spac- 
ings. It turns out that cq value is practically independent 
upon the concrete choice of the fitting range being 1.03(2) 
and 1.05(2) in both cases correspondingly. Therefore the 
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value of the string tensions ratio in the continuum limit 
is 



jSU{2) 



1.04(3). 



(37) 



Note that this result is indeed remarkable. Taken at face 
value it indicates that in the continuum limit the full 
SU{2) string tension arises solely due to the rather weak 
HP^ projected fields. Simultaneously this is the first in- 
dication that what had been cut away from the original 
fields by HP^ projection corresponds likely to pure per- 
turbation theory. On the other hand, Eq. (|37l) should not 
come completely unexpected. Indeed, the non-vanishing 
value of the dimension 2 condensate in HP^ projected 
fields (see section IV A|l is a hint that the projected the- 
ory might be confining 0, 0, ■ However, we will not 
dwell on this issue any longer but consider the anatomy 
of HP^ string tension instead. 



2. String Tension and Topological Fluctuations 

One of the central point of Ref. Jal was rediscovery 
of the lumpy structure 143. lisl l4a[ of the topologi- 
cal charge density bulk distribution in HP^ cr-model ap- 
proach. As was stressed in that paper, the term 'lumps' 
is in fact uncertain until one introduces a particular cut- 
off Kq on the magnitude of the local topological charge 
density. Indeed, the most straightforward argument here 
is that in the numerical simulations the topological den- 
sity is always known with finite accuracy and therefore a 
particular cutoff applies inevitably. Moreover, the intro- 
duction of finite is inherent to practically all studies of 
the gauge fields topology. For instance, the overlap-based 
topological density, which is given by the sum of Dirac 
eigenmodes V'a contributions, is usually either restricted 
to lowest modes, A < A, or is considered for all modes 
available on the lattice, A < 1/a. In both cases it is ap- 
parent that a particular cut on the topological density is 
introduced although it might not be simply expressible 
in terms of A or 1/a. Another example of this kind is 
provided by the investi gati ons of local chirality of Dirac 
eigenmodes li^ l43.T45ll4^ . Here the local chirality 
is considered only in points at which the eigenmode is 
reasonably large. 

Then the crucial question is the dependence of phys- 
ical observables upon Kq. The point of view accepted 
in Ref. '3| is that this dependence must be trivial since 
Kq is technical rather than physical parameter. This 
requirement fixes, in fact, the physically meaningful 
values of Kq ('physical window' in the terminology of 
Ref. 01) which appears to be rather narrow and is around 
[300 - 350 MeV]"*. The attempt to consider the limit 
Ag ^ results in rather peculiar structure of topological 
excitations which is reminiscent to one discovered in [J 
with overlap-based definition of the topological density. 
Namely, the density of lumps seems to be divergent with 



Kq ^ Q due to the abundance of the lumps consisting 
of only one lattice site. Simultaneously lumps of larger 
volumes are accumulating in a few (typically two) perco- 
lating lumps with linear extent equal to the size of the 
lattice. However, both the small lumps distribution and 
the volume occupied by percolating lumps are strongly 
Aq-dependent and for that reason the small values of the 
cutoff were excluded in Ref. from the Ag physical win- 
dow. It is important that the topological susceptibility 
remains lumps-saturated even beyond the above perco- 
lation transition and essentially this was the justification 
for the existence of Ag physical window. 

It is natural then to ask the same question about the 
HP^ string tension, namely, to consider Ag dependence 

of cr^P , Eq. 1(23). Evidently we should confront the 
Ag physical window found in with the requirement 

of Ag independence of a™' . The first problem to be ad- 
dressed here is how the string tension could at all be 
dependent upon the cutoff on the topological density. 
For the topological susceptibility the solution is trivial, 
it suffice to consider the small values of the topologi- 
cal density \q{x)\ < Kq as being zero exactly. Now we 
have to formulate the same principle at the level on HP^ 
projected gauge fields maintaining the gauge covariance. 
The straightforward way to achieve this is provided by 
Eq. (O. Indeed, every HP^ induced link is representable 
as 



U, 



HP^ 



\{qx\qx+^L}\ 



(38) 



\1 + LO^LO^+f^l 



'!Jx-\-fi — Sx ^x.^ 9x 



+1-1 J 



where h^^^, € SU{2) is, in fact, gauge invariant quan- 
tity. If the topological density in the vicinity of point 
X is small, \q{x)\ < Kq, it is sufficient to equate the in- 
homogeneous coordinates tUx, Wx+n to each other thus 
replacing 



1 



C/, 



gx 9x^ 



(39) 



This receipt is indeed natural since it maintains gauge 
covariance and nullifies the contribution of the link U^^^ 
to any Wilson loop. Note however that we are not al- 
lowed to change the links and topological density at other 
points. This means that Eq. should be understood 
on the level of links only, it is not allowed to actually 
change neither | q^ ) nor | qx+n ) ■ Our implementation of 
this approach is as follows: the link matrix is trans- 

formed according H39|) if the magnitude of the topological 
density in 8 neighboring hypercubes containing the link 
is smaller than Ag. 
Unfortunately, the only set of configurations on which 
we could analyze Ag dependence of is the one cor- 
responding to /3 = 2.475 on 16'' lattice, thus no scahng 
checks are possible at the time being. The Ag dependence 

of the ratios a^^\Kq) / a^^' {Q) and v/x(Ag)/x(0), 
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FIG. 11: Ag dependence of the HP^ string tension (squares) 
and of the topological susceptibility (circles) both normalized 
to their respective values at Ag = {(3 — 2.475, 16*). Shaded 
box stands for the approximate location of the lumps perco- 
lation transition (see text). 



lattice) regions of sign-coherent topological charge first 
discovered in and discussed within HP^ approach 
in 5j. Therefore the string tension in HP^ projected 
configurations, which accounts for the full SU{2) string 
tension in the continuum limit, is most sensitive not to 
the magnitude of the topological density fluctuations, but 
rather to the long range correlations between them. Con- 
trary to that the topological susceptibility is almost blind 
to such correlations and, in fact, is saturated by almost 
random lumps in the topological density bulk distribu- 
tion. 



3. Polyakov lines correlation function 

Let us discuss another conventional confinement in- 
dicating observable, namely, the Polyakov line and the 
corresponding correlation function. The basic observa- 
tion here is that the Polyakov line ceases to be an order 
parameter once the gauge potentials (or link matrices in 
lattice settings) are written as 



where cr^^^ (0) and x(0) are the HP^ string tension and 
topological susceptibility with no cutoff imposed, is pre- 
sented on Fig. ^2 It is remarkable that the sensitivity 
of these quantities to the Ag cutoff are radically distinct. 
The susceptibility decreases only slightly with rising 
and is changed by « 20% in the whole range of cutoff con- 
sidered. Contrary to that the HP^ string tension drops 
down to zero at A^-a^ « 3-10^'* and remains vanishing for 
all higher values of the cutoff. Moreover, we cannot ex- 
clude the possibility of the phase transition in the infinite 
volume limit. It is crucial that the string tension vanishes 
just below the Ag physical window Kq-a^ G (4-;- 8) • 10"'* 
discussed in Ref. 5]. In turn this is not surprising since 
the HP* string tension was not considered in that paper. 

In fact, it is not difficult to identify what happens phys- 
ically around the point Ag • a** « 3 • lO"**. Indeed, exactly 
around this value of Ag the percolating lumps cease to 
exist, the gap separating percolating lumps from the re- 
maining spectrum in the lumps volume distribution dis- 
appears. To illustrate this point we indicated the ap- 
proximate location of the lumps percolation transition 
by shaded box on Fig. Unfortunately, at present we 
are not able formulate this qualitative picture more rig- 
orously, mainly because of the limited data we have right 
now. Nevertheless, we do believe that the disappearance 
of HP* string tension could be formulated equivalently in 
the language of lumps percolation transition. The cor- 
responding strong correlation is evident already on our 
limited data sets and it is probably the question of rel- 
atively short time to make this conjecture numerically 
convincing [l6j . 

To summarize, the investigation of HP* string ten- 
sion dependence on the cutoff imposed on the topologi- 
cal charge density reveals that the string tension is most 
likely due to the prolongated (extending through all the 



( Qx I gx+M ) 

\{qx\qx+t,)\ 



(40) 



The intrinsic reason for this is that the relevant global 
center symmetry cannot be formulated in terms of | (7^; ) 
fields. Hence the expectation value (P) of the Polyakov 
line is not obliged to vanish in confinement phase once the 
representation (|4()|l is adopted. Physically this is because 
it is always possible to form the gauge invariant compos- 
ite state \qx)i^x from \q) field and quark operator -0. 
More explicitly, consider for instance the derivation pre- 
sented Ref. Wfl^. The equation which is crucial here is 
the static time-evolution equation 



{dt + Ao)^{t,x) = 0, 



(41) 



which eventually produces the Polyakov line gauge trans- 
porter to be averaged in the thermal ensemble (note that 
our gauge potentials are anti-hermitian, hence there is no 
imaginary unit in (|41|l ). Within the representation l|40l) 
this equation is equivalent to (we omit the inessential x 
dependence here) 



[qt\dt[\qt)^t 



0, 



(42) 



which could be solved as 
t 

^t = {qt\ j[\q^)-f] + (gtko)-V'o, (43) 



where the state | (7^ ) is orthogonal to \q), {q\q^) = 0, 
and is determined, in fact, uniquely, / is arbitrary func- 
tion and the second term on the r.h.s. is determined from 
boundary condition at t = 0. Evidently, for single quark 
free energy Fq, considered in the canonical formalism. 



14 



only the last term does matter (see Ref. '4^ for details) 
and we obtain 



= Tr[e-«/^ 



(44) 



(9i/Tko)] = 1, 



where H is the Hamiltonian, T is the temperature and we 
assumed that the periodic boundary condition in imagi- 
nary time are imposed on the | q ) fields. The result is in 
agreement with what was noted above: the color charge 
of both ijji and t/iq gets screened by the quanta of | g ) 
fields. Note that Eq. (|44() is not dynamical, it is valid 
irrespectively at both low and high temperatures. More- 
over, it confirms that the Polyakov line expectation value 
(P) measured on HP^ projected fields is not related to 
the free energy of static color source and generically is 
not vanishing at low temperatures 



(P) ^ 0. 



(45) 



In fact, Eq. 1)45(1 is strictly confirmed by our measure- 
ments. Depending on the lattice size and spacing (P) 
varies from « 0.42 to « 0.67 revealing non-trivial depen- 
dence on both V and a. However, we did not investigated 
this dependence due to the uncertainty of its physical im- 
plications. On the other hand, the physical meaning of 
the Polyakov lines connected correlation function is not 
changing 



(P) 



1 



cxp{--^"^ (|x|)}. 



(46) 



Due to the weakness of HP^ projected fields the behavior 
(I46|l is firmly confirmed even on the largest lattices we 
have, ft turns out that the potential V^^ (R) extracted 
from (|46|l is in agreement with what we have discussed 
in section IV B II although it has much larger error bars 
due to the subtraction of disconnected contribution. In 
particular, it approaches linear asymptotic at large dis- 
tances, the corresponding string tension being in agree- 
ment with that of section IVB II 



C. Chiral Fermions in HP^ Projected Gauge 
Background 

It had long been understood that the problem of con- 
finement and chiral symmetry breaking are closely re- 
lated 48, 49] (for recent discussions see, e.g., Refs. |50|). 
In this section we investigate the properties of HP'^ pro- 
jected fields (|30|) as they are seen by fermionic probes 
provided by the overlap Dirac operator (|23|) . Physically 
the most important part of the spectrum of ((23(1 is given 
by low lying eigenmodes 



(47) 



the spectral density p(A) of which gives the quenched 
chiral condensate via Banks-Casher relation IsJ 
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FIG. 12: The spectral density of lowest Dirac eigenmodes 
calculated on fixed volume configurations, Table|I] at diflferent 
spacings. Line represents linear interpolation towards A = 0. 



('0'0) — liniA^o p(A). Another problem of prime im- 
portance is the localization properties of low lying modes, 
which is conventionally formulated in terms of the inverse 
participation ratio (IPR) I\ (see, e.g., Refs. JJ3,, 11, sj] 
for definitions). Specifically, we are interested in the scal- 
ing of fermionic IPRs with lattice spacing and will not 
consider dependence upon the total volume (see Refs. [s^ 
for the similar treatment of scalar probe particles). We 
study both p(A) and Ix on HP^ projected gauge fields 
obtained from our fixed volume data sets (Tabled which 
cover rather wide range of lattice spacings and allow to 
investigate the scaling of p(A) and I\ in details. Note 
that the spacing dependence of IPR for near zero eigen- 
modes is still the subject of intensive discussions in the 
literature TO, 11, 52]. Thus the HP^ cr-model embedding 
method again provides the unique opportunity to study 
the issue. 

One technical note is now in order. To speed up the 
calculations we considered only 20 lowest Dirac eigen- 
modes and excluded the exact zero modes from the anal- 
ysis since zero modes are known to be physically quite 
different from the remaining part of the spectrum. This 
means that we were able to collect the sufficient statistics 
only for A < 400 MeV and for this reason all the graphs 
below are snipped off for higher eigenvalues. However, 
this does not spoil the significance of our results since 
exactly the region A < 400 MeV is most important phys- 
ically. 

Let us consider first the spacing dependence of the 
spectral density px at small A as is measured on HP^ 
projected gauge fields. Fig. E| Note that for readabil- 
ity reasons not all available graphs are depicted since the 
data sets being expressed in physical units look essen- 
tially the same at different spacings. It is clear that the 
points at various a are falling, in fact, on the top of each 
other thereby confirming that the spectral density of low 
lying modes is spacing independent. This not only jus- 
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FIG. 13: IPRs Ix for lowest eigenmodes at various lattice 
spacings. Note the explicitness of the 'mobility edge' around 
A ^ 200 MeV. 




FIG. 14: IPR for lowest eigenmodes, Eq. I|49|l . as a function 
of lattice spacing. The fitting curves 15011 are represented by 
various lines (see text). 



tifies the assertion made earlier that the HP^ projection 
captures correctly the topological aspects of the original 
gauge fields, but also allows to estimate the quenched 
chiral condensate 



iih'it') = TT lim piX) 



= [278(6) MeV]^ 



(48) 



which is in agreement with its value measured on the 
same original gauge configurations in Ref. IQj and is 
fairly compatible with magnitude of (■(/'V') known in the 
literature (see, e.g., Refs. w| and references therein). 

As far as the localization properties of lowest eigen- 
modes are concerned , ou r results a re q ualitatively similar 
to ones reported in (see also ^^H^)- Namely, for 
A < 200 MeV the eigenf unctions show notable localiza- 
tion which is apparent from Fig. 1131 the corresponding 
IPRs are ranging from 4-=-8 at A w 150 MeV to 10-^25 at 



smallest non-zero eigenvalues. However, the most impor- 
tant problem is the spacing dependence of I\, to study 
which it is advantageous to consider the IPR averaged 
over all lowest modes < A < A = 50 MeV 



1 



(49) 



A<A 



where N/^ is the total number of modes in the above 
interval. This allows to notably improve the statistics, 
however, we checked that the scaling properties of / are 
A-independent provided that A < 100 MeV. The lattice 
spacing dependence of / allows to investigate quantita- 
tively the dimensionality d of the regions which localize 
the lowest eigenmodes. As it is customary now (see, e.g., 
Refs. [13 111 El for details) we fitted the I{a) depen- 
dence to 



/(a) = bo + hi 



(50) 



the best fitting curves for = 0, 1, 2, 3 are shown by lines 
on Fig. El Unfortunately, our data does not allow to 
pinpoint the dimensionality d with high enough accuracy. 
However, it still favors the two-dimensional structure of 
the localization regions although we cannot fairly exclude 
d = 1 possibility. The relevant /n.d.f. is equal to 0.33 
for d — 2 and 0.62 for d — \, while being significantly 
larger in c? = and d = 3 cases (1.26 and 0.82 corre- 
spondingly). In either case the data is incompatible with 
naively expected localization on four-dimensional 'balls' 
of finite physical extent, the localization volumes are def- 
initely shrinking in the limit a — > 0, which is in accord 
with what is known in the literature [lol [Til Is^ about 
the low lying fermionic eigenmodes in equilibrium gauge 
background. 



D. Summary 

Let us briefly summarize the results obtained so far on 
the dynamics of HP^ projected gauge fields. We stress 
from very beginning that the HP^ projection is com- 
pletely different from what is usually known by the term 
'projection' in the literature. Namely, it maintains the 
gauge invariance and in view of the investigated irrele- 
vance of Gribov copies problem provides unique gauge in- 
variant variables characterizing the original gauge back- 
ground. The main conjecture of this section to be jus- 
tified below is that HP"^ projected fields represent al- 
most entirely the non-perturbative content of the orig- 
inal gauge fields. We illustrated this point from various 
perspectives starting from the simplest possible observ- 
able, namely, the HP^ induced gauge curvature. Remark- 
ably enough, the projected fields appear to be extremely 
weak, the averag ed plaquette (1/2 TrJ/^^ ) being dis- 
tributed with striking power law, which is to be com- 
pared with usual exponential distribution. As a conse- 
quence, the lattice spacing dependence of (1/2 Trf/^^^ ) 
is totally different from what could be usually expected. 
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we found no sign whatsoever of the perturbation theory 
in it. As a matter of fact 1 — (1/2 Trf/^^ ) vanishes in 
the hmit a ^ 0, the leading spacing corrections being 
of order 0(a'*) and O(a^) only. The corresponding co- 
efficients are to be identified with the gluon condensate 
and quadratic power correction, both of which we mea- 
sured with highest accuracy available so far. As far as 
the gluon condensate is concerned, its value is in a com- 
plete agreement with the existing literature. At the same 
time, the numerically precise value of the quadratic cor- 
rection given in lattice spacing normalization is reported 
for the first time and allows to pose the problem of un- 
usual power corrections in Yang-Mills theory on qualita- 
tively new level. 

Then we addressed the problem of confinement in the 
HP^ projected gauge fields. Despite of the relative weak- 
ness of the HP^ induced potentials we showed that the 
projected theory is still confining. Although at finite lat- 
tice spacing the projected string tension is smaller than 
that on the original fields, it nevertheless reproduces the 
full SU{2) string tension in the limit a ^ 0. More- 
over, the last statement was found to be quite robust 
with respect to the technicalities involved. However, it 
would be plainly wrong to rely on the HP^ projection at 
all distances in particular because the method is prac- 
tically blind to the perturbation theory. It turns out 
that the heavy quark potential at small R < 0.2 fm 
scales clearly violates the reflection positivity require- 
ments, which, however, is to be expected due to the 
intrinsic non- locality of HP^ cr-model embedding. We 
also argued that within our approach the usual confine- 
ment sensitive observables like Polyakov line expectation 
value (P) is to be considered with care. In particular, 
the notion of global center symmetry, which is crucial to 
ensure (P) = at low temperatures in usual settings, 
does not exist within the HP^ projection method. This 
clearly rises the question on the relevance of the center 
symmetry to confinement , ho wever, we postponed the 
corresponding discussions |l6j | and concentrated only on 
the experimental results instead. 

Since the HP^ projection is ultimately related to the 
topology of the gauge fields, it makes possible the inves- 
tigation of the microscopic origin of the projected string 
tension. We found that the global (percolating) regions 
of sign-coherent topological charge 0, 0| seem to be fun- 
damentally important. Namely, the geometrical explicit- 
ness of HP^ cr-model, being superior to any other known 
topology considerations, allows to scan the gauge fields 
topology at various scales and to associate rather unam- 
biguously the projected string tension with the existence 
of the above regions. Note that at present much more 
data is needed to put the issue on the numerically solid 
grounds, but for us the correlation is too much evident 
to be in any doubts. 

Finally, in order to present the comprehensive picture 
of the projected fields dynamics we considered the spec- 
trum of the overlap Dirac operator in the HP^ induced 
background. The observables here are the spectral den- 
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sity of low lying eigenmodes and their localization prop- 
erties conventionally encoded into inverse participation 
ratio (IPR). We found that both this quantities are es- 
sentially the same as they were on unprojected fields al- 
though the actual measurements are much simpler due 
to the weakness of HP^ potentials. 

However, the above conclusions rely strongly on the 
fact that HP^ projection leaves aside only the pertur- 
bation theory. It is the purpose of the next section to 
demonstrate that indeed what had been left after the 
HP^ projection likely corresponds to purely perturbative 
fields. 



VI. WHAT HAD BEEN CUT AWAY BY HP^ 
PROJECTION? 

The problem to be considered below is the properties 
of the fields Yx.ii which are the part of original configu- 
rations cut away by HP^ projection (|30|1 . Symbolically 
Vx.iL could be defined by = U jU^^ , however, this 
equation makes no sense as it stands. Indeed, both V 
and J7^^ transform as usual under the gauge transfor- 
mations and hence the transformation properties of Y are 
undefined. It is clear that we could make sense of Y only 
in a particular gauge, however, we are not aware of any 
detailed treatment of this generic problem in the lattice 
context. The only exception is probably Ref. [ISj where 
related issues are discusses. Although the extraction of 
Y variables seems to be similar to the background fields 
technique, we note that usually the background poten- 
tials are taken for granted already in a particular gauge 
thus invalidating any similarities with the present case. 

In order to provide the unambiguous meaning to T4^^ 
fields we fixed the most convenient and well defined 
gauge, in which all | qx ) variables are given by the r.h.s. 
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of Eq. © 



Then Vx^^ variables are defined by 



• [ 1 ; 



Ux.f_L 



u. 



HP^ 



(51) 



(52) 



Apparently it is important here that U^^^ is expressed in 
terms of factually gauge invariant inhomogeneous coordi- 
nates Wa:, ujx+^,- Note that the results to be discussed be- 
low are dependent in principle upon the particular gauge 
chosen, at least we don't have precise theoretical argu- 
ments to exclude this possibility. Since we agreed from 
very beginning not to go deep into theoretical analysis 
we leave this problem until the future work 16]. 

In this section we consider the topological aspects of 
V fields (|52|l . investigate their confining properties and 
the spectrum of overlap Dirac operator in V fields back- 
ground. For this purposes we generated 520 configura- 
tions of V fields from the original data sets at /3 = 2.400 
and (3 ~ 2.475 on 16^ lattices. However, let us start 
again from the simplest observable 1/2 TrV^ and con- 
sider its distribution compared to that of the original 
fields. The typical graph of this sort, taken on a particu- 
lar (3 — 2.475, 16^^ configuration, is presented on Fig.lTSl 
In fact, for all our V configurations the distributions of 
l/2TrVp look very similar. Remarkably enough, it is 
hardly possible to distinguish the points belonging to the 
original gauge potentials and to V variables, the only dif- 
ference is seen very close to unity and results in < 1% 
larger value of (1/2 Tr Vp) on this configuration compared 
to that of (l/2Tr?7p). However, it should not be sur- 
prising in view of the above discussed weakness of HP^ 
projected fields. It might be tempting to conclude that 
the V fields are essentially the same as the original gauge 
potentials and then the validity of HP^ projection would 
become problematic. However, it is the purpose of this 
section to show that V fields are completely trivial as far 
as the non-perturbative aspects are concerned. 



A. Triviality of Topological Properties 

The technically simplest test of topology is provided 
by the application of HP^ cr-model embedding method 
to V configurations. For us it was almost shocking to 
discover that the Qhr^ topological charge is identically 
zero for all V configurations we have generated on var- 
ious lattices. Of course, Eq. H52I) expresses our specific 
intent to cut away the topology related aspects of the 
original configurations, however, the identically vanish- 
ing topological charge is indeed remarkable. Evidently 
this must be checked by other methods to ensure that 
Q = is not the artifact of HP^ approach. We did 
this cross-check with overlap-based definition and indeed 




FIG. 16: — ln{W(r, R)) as a function of T at various R mea- 
sured on V configurations generated on /3 = 2.475, 16* lattice. 



found that Q overlap is equal to ±1 only on two V con- 
figurations out of 520. Evidently, this disagreement of 
Qnpi and Qoveriap is completely inessential and should 
be attributed to what had been discussed in section ITTll 
Therefore we are confident that the topological proper- 
ties of V fields are entirely trivial, the topological charge 
and susceptibility vanish identically 



Q = 0, 



X 



0. 



(53) 



This is the first indication that what had been left aside 
by HP^ projection is the pure perturbation theory. 



B. Loss of Confinement 

The next question to be addressed is the possibility to 
have confinement on V configurations despite of the triv- 
ial topology inherent to them. Note that the answer is 
by no means evident since the relation between topology 
and confinement is still not proved rigorously. Consider 
first the Wilson loops confinement criterium. We mea- 
sured the expectation values of planar TxR Wilson loops 
on our V configurations in exactly the same way as was 
done in section IVB II As a matter of fact, on all our V 
data sets the corresponding heavy quark potential fiat- 
tens at large distances and the string tension vanishes. 
In order to convince the reader that the potential is in- 
deed constant at large separations we plotted on Fig. ^| 
the expectation values — \n{W{T, R)) at various i? as a 
functions of T. It is apparent that the slope of the curves 
remains constant for i? > 3 although it is still rising at 
smaller distances. The corresponding potential (Fig. I17() 
was fitted to analytical equation analogous to l|35|l and 
turned out to be totally compatible with Coulomb law 



ViR) 



const 



a 
R 



(54) 
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FIG. 17: The heavy quark potential extracted from V fields 
on /3 = 2.475, 16"* lattice. The curve is the fit to analytical 
dependence (1541 . 



FIG. 18: The spectral density of overlap Dirac operator mea- 
sured on V fields. Note the complete absence of eigenmodes 
in A < 200 MeV interval. 



However, we do not presenting the best fitted value of 
a parameter since on our lattices it would be dictated 
mostly by V{R) at smallest R, which is highly sensitive 
to the details of smearing and hypercubic blocking used 
in our calculations (note that essentially for this reason 
the points R = 1,2 deviate from 154|) '). 

As far as the Polyakov line and the corresponding cor- 
relation function are concerned, the arguments of sec- 
tion IVB 31 are clearly invalid for V fields and therefore 
both (P) and (PqPx) are expected to be the confinement 
sensitive observables. It turns out that the expectation 
value (P) is strictly positive on all our V data sets be- 
ing essentially the same at /3 = 2.400, (P) = 0.294(2), 
and at /3 = 2.475, (P) = 0.290(4). The strict positivity 
of Polyakov line distributions apparently rises the ques- 
tion on the relevance of center symmetry, however, due 
to the possible lack of statistics we'll not dwell on this 
issue. The Polyakov lines connected correlation function 
{PqPx) — {Py measured on V configurations also clearly 
shows deconfining behavior, namely, it falls off only as a 
power of I a; I being strictly linear on log-log plot. Note 
however that the subtraction of the disconnected part 
violently spoils the available numerical accuracy and for 
this reason we do not quote the results of the fits. 

To summarize, we found clean evidences that what 
had been cut away by HP^ projection corresponds to 
deconfining theory, in which the heavy quark potential 
is totally compatible with Coulomb law. This provides 
the most stringent illustration that the V fields configu- 
rations correspond to pure perturbation theory with no 
sign whatsoever of neither non-trivial topology nor con- 
finement. 



C. Fermionic Probes in V Background 

To conclude this section let us consider the properties 
of low lying eigenmodes of overlap Dirac operator in V 
fields background. The measurements performed here are 
essentially the same as was done in section IV CI Note 
that we considered only V data sets corresponding to 
/3 — 2.475 for reasons to become clear shortly. 

The spectral density of low lying eigenmodes is pre- 
sented on Fig. [TBI and it is apparent that the spectrum is 
completely empty for A < 200 MeV, no one eigenmode 
was found in this region on all available V configurations 
at /? = 2.475. Essentially this precludes any considera- 
tions of low lying eigenmodes, there are strictly speaking 
no modes to be discussed here. Thus the spectral den- 
sity for small A is identically zero leading to vanishing 
quenched chiral condensate 

i^pip) = 0. (55) 

Evidently it makes no sense to discuss the localization 
properties of low modes in V background, the quantities 
like I\ or /, Eq. H49|l . become undefined. Nevertheless, it 
is still instructive to qualitatively consider the change in 
the degree of localization for eigenmodes at A > 200 MeV 
compared to that on the original and HP^ projected fields 
(see section IVCII . It turns out that the degree of local- 
ization expressed by IPRs I\ drop down significantly and 
for all calculated modes A > 200 MeV it is around 1.2, 
which essentially means complete delocalization. 



D. Summary 

In this section we presented rather clean and concise 
evidences that what had been left aside by HP^ projec- 
tion H30|l corresponds to purely perturbation theory. The 
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primary aim of this presentation was to show that the re- 
sults of section |V] are not invalidated by the components 
Vx,fi of the original gauge fields which do not survive 
HP^ projection. We considered this issue from various 
perspectives and showed that V fields are not confining, 
the corresponding heavy quark potential is adequately 
described by usual Coulomb law. Moreover, the topolog- 
ical properties of V configurations are completely trivial, 
leading, in particular, to the identically vanishing topo- 
logical susceptibility. This result turns out to be quite 
robust with respect to different definitions of lattice topo- 
logical charge. As a consequence the low lying spectrum 
of overlap Dirac operator in V background has nothing in 
common with that on the original or HP^ projected fields. 
Strictly speaking, the low lying spectrum does not exist, 
no one eigenmode was found below A < 200 MeV despite 
of rather large number of analyzed configurations. This 
means that quenched chiral condensate vanishes identi- 
cally in V background, {'tpip) = 0. 

Therefore, we are confident that the content of the 
original fields, which is cut away by HP^ projection, is 
completely trivial as far as the non-perturbative aspects 
are concerned. This provides at least numerically solid 
grounds for the results obtained in section IVI 



VII. CONCLUSIONS 

The primary goal of the present publication was the 
further development of HP^ a-model embedding ap- 
proach aimed to investigate the topology of SU{2) 
gauge fields. While Ref. discussed in details the the- 
oretical aspects, the tests of actual numerical implemen- 
tation were not so much convincing. In this paper we 
filled this gap and presented the high statistics com- 
parison of HP^ based topological charge, QhpIj with 
ones obtained via field-theoretical, Qdover, and overlap 
Dirac operator based, Qoveriap, approaches. As far as the 
field-theoretical definition is concerned, the correspond- 
ing comparison is meaningful only on semiclassical con- 
figurations and in this case we essentially found the com- 
plete equivalence of Qnpi and Qdover- In the most inter- 
esting case of equilibrium vacuum configurations we con- 
fronted Qnpi and Qoveriap On rather large number of data 
sets at different spacings and found that both definitions 
are strongly correlated and identify essentially the same 
topology. Unfortunately, it was difficult to compare our 
results with similar calculations in the literature, how- 
ever, the correlation between Qupi and Qoveriap turns 
out to be quantitatively the same as it was in the case 
of Qoveriap and Qdover definitions. Basing on that we 
concluded that the HP^ cr-model embedding method is, 
in fact, computationally superior to any other definitions 
of the lattice topological charge available so far. 

The next immediate problem is to investigate on the 
numerically confident level the scaling of the topologi- 
cal susceptibility both with lattice spacing and volume. 
Besides of the methodological importance, this problem 



is crucial in order to address the physical relevance of 
the method itself and its sensitivity to lattice artifacts. 
We performed statistically convincing scaling checks of 
the topological susceptibility using the HP^ embedding 
approach and found almost perfect scaling with both lat- 
tice spacing and volume. The topological susceptibility 
being extrapolated to the continuum limit turns out to be 
^^^/'^ = 216(4) MeV and is totally compatible with what 
is known in the literature. The scaling properties of x 
make us confident that the HP^ embedding method is 
insensitive to lattice dislocations, moreover, the topolog- 
ical susceptibility does not require, in fact, neither mul- 
tiplicative nor additive renormalizations. 

As might be anticipated from the results of Ref. 0, 
the HP^ embedding approach is not restricted to the 
topology related measurements. Its essence is to pro- 
vide a unique (modulo Gribov copies problem which, 
however, was proved to be inessential) configuration of 
quaternionic valued scalar fields to any given original 
gauge background. Then the corresponding HP^ projec- 
tion comes out naturally and we considered in details the 
properties of HP^ projected fields. Note that HP^ projec- 
tion is radically different from what is commonly known 
by this term, in particular because it maintains the gauge 
invariance from very beginning. The first striking obser- 
vation here is that the mean curvature of projected po- 
tentials is extremely small, while its spacing dependence 
shows no sign whatsoever of the perturbation theory. In 
fact, it only contains the terms which are to be iden- 
tified with giuon condensate and non-trivial quadratic 
correction to {asG^/n). We calculated both this quan- 
tities with inaccessible so far accuracy, the result for 
the gluon condensate being {as G^/tt) — 0.066(2) GeV* 
which fairly agrees with the existing literature. On the 
other hand, the numerically precise value of the quadratic 
correction satisfies the known tight bounds and allows to 
pose the problem of unusual power corrections in Yang- 
Mills theory on qualitatively new level. 

Then we addressed the crucial problem of confinement 
in the language of HP^ projection. It is remarkable that 
despite of the weakness of the projected fields the HP^ 
string tension is not only non-zero at finite lattice spac- 
ing, it accounts, in fact, for full SU(2) string tension in 
the continuum limit, \ima-,o[(J^^' /a^'^^'^'^]^/'^ = 1.04(3). 
Moreover, the explicitness of the gauge fields topology in 
HP^ variables allows us to study the microscopic origin 
of the projected string tension. We argued that ^ 
is to be associated with the existence of percolating re- 
gions of sign-coherent topological charge, first discovered 
(with overlap-based approach) in and then confirmed 
in HP^ method in As far as we can see, this is the 
first explicit demonstration of the confinement-topology 
interrelation. 

We concluded the investigation of HP^ projected fields 
with consideration of low lying spectrum of overlap Dirac 
operator in the projected background. The spectral den- 
sity of low eigenmodes turns out to be essentially the 
same as it was on the original gauge fields and allows us 
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to estimate the value of chiral condensate in the quenched 
approximation {■ip^)^/'^ — 278(6) MeV, which is fairly 
compatible with its commonly accepted value. Simulta- 
neously we measured the localization properties of low 
eigenmodes, which are conventionally expressed in term 
of the inverse participation ratio (IPR). Knowledge of 
IPRs at various spacings allows us to conclude that the 
localization regions are likely to be either two- or three- 
dimensional. 

It was understood from very beginning that the valid- 
ity of the above results relies strongly on the ability of 
HP^ projection to capture only non-perturbative aspects 
of the original gauge fields. Therefore it is mandatory 
to study what had been left aside by HP^ projection. 
We investigated carefully this issue by showing first that 
the part of the gauge fields not surviving HP^ projection 
could be defined unambiguously. Then it is the matter 
of straightforward calculation to show that indeed what 
had been cut away by HP^ projection is completely triv- 
ial as far as the non-perturbative aspects are concerned. 
In particular, there is no topology in these fields, the 
topological charge and susceptibility vanish identically; 
confinement is definitely lost although the Coulomb part 
of the potential is clearly seen; and finally the low Dirac 



spectrum completely disappears, not a single low eigen- 
mode was found on our rather large set of configurations. 
Thus we concluded that the HP^ projection indeed cap- 
tures only the non-perturbative content of gauge back- 
ground and what is left aside corresponds to pure per- 
turbation theory. 

Finally let us remind our original intent not to discuss 
the theoretical implications of the above findings. Only 
experimental numerical results had been reported since 
otherwise the discussions would definitely lead us far be- 
yond the scope of present publication. The corresponding 
theoretical treatment is to be published elsewhere p^ . 
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